On the effect of nocturnal radiation fog on the development of the daytime convective boundary layer: A large-eddy simulation study

The potential effect of failing to predict nocturnal deep radiation fog on the development of the daytime convective boundary layer (CBL) is studied using large-eddy simulations. Typical spring and autumn conditions for the mid-latitudes are used to perform simulations in pairs. Fog formation is allowed in one simulation of each pair (nocturnal fog [NF]) and is suppressed in the other (clear sky [CS]). This allows for the identification of properties (temperature, humidity, boundary-layer depth), conditions, and processes in CBL development that are affected by fog. Mixing-layer temperatures and boundary-layer depths immediately after fog dissipation in CSs are shown to be up to 2.5 K warmer and 200 m higher, respectively, than the NF counterparts. Additionally, greater water vapor mixing ratios are found in the CSs. However, owing to greater temperatures, relative humidities at the CBL top are found to be less in CSs than in the corresponding NFs. This relative humidity difference might be an indication that cloud formation is suppressed to some extent. The magnitude of the differences between CSs and NFs during the day is mainly correlated to the fog depth (in terms of duration and liquid water path), whereas the key processes responsible for differences are the atmospheric long-wave cooling of the fog layer (for temperature development) and droplet deposition (for water vapor mixing ratio development).publishedVersio

An Investigation of the Grid Sensitivity in Large-Eddy Simulations of the Stable Boundary Layer

We revisit the longstanding problem of grid sensitivity, i.e., the lack of grid convergence in large-eddy simulations (LES) of the stable boundary layer. We use a comprehensive set of LES of the well-known Global Energy and Water Cycle Experiment Atmospheric Boundary Layer Study 1 (GABLS1) case with varying grid spacings between 12.5 m and 1 m to investigate several physical processes and numerical features that are possible causes of grid sensitivity. Our results demonstrate that there are two resolution regimes in which grid sensitivity manifests differently. We find that changing the numerical advection schemes and the subgrid-scale models alters the simulation results, but the options tested do not fully address the grid-sensitivity issue. Moreover, sensitivity runs suggest that the surface boundary condition and the interaction of the surface with the near-surface flow, as well as the mixing with the free atmosphere, are unlikely to be the causes of the observed grid sensitivity. One interesting finding is that the grid sensitivity in the fine grid-spacing regime (grid spacings ≤2m) is closely related to the reduction in the energy content of large-scale turbulence, leading to less turbulence kinetic energy and hence lower boundary-layer heights. The present work demonstrates that there is still an urgent need to address this grid-sensitivity issue in order to perform reliable LES of the stable boundary layer.publishedVersio

Towards a better representation of fog microphysics in large-eddy simulations based on an embedded Lagrangian cloud model

The development of radiation fog is influenced by multiple physical processes such as radiative cooling and heating, turbulent mixing, and microphysics, which interact on different spatial and temporal scales with one another. Once a fog layer has formed, the number of fog droplets and their size distribution have a particularly large impact on the development of the fog layer due to their feedback on gravitational settling and radiative cooling at the fog top, which are key processes for fog. However, most models do not represent microphysical processes explicitly, or parameterize them rather crudely. In this study we simulate a deep radiation fog case with a coupled large-eddy simulation (LES)–Lagrangian cloud model (LCM) approach for the first time. By simulating several hundred million fog droplets as Lagrangian particles explicitly (using the so-called superdroplet approach), we include a size-resolved diffusional growth including Köhler theory and gravitational sedimentation representation. The results are compared against simulations using a state of the art bulk microphysics model (BCM). We simulate two different aerosol backgrounds (pristine and polluted) with each microphysics scheme. The simulations show that both schemes generally capture the key features of the deep fog event, but also that there are significant differences: the drop size distribution produced by the LCM is broader during the formation and dissipation phase than in the BCM. The LCM simulations suggest that its spectral shape, which is fixed in BCMs, exhibits distinct changes during the fog life cycle, which cannot be taken into account in BCMs. The picture of the overall fog droplet number concentration is twofold: For both aerosol environments, the LCM shows lower concentrations of larger fog droplets, while we observe a higher number of small droplets and swollen aerosols reducing the visibility earlier than in the BCM. As a result of the different model formulation we observe higher sedimentation rates and lower liquid water paths for the LCM. The present work demonstrates that it is possible to simulate fog with the computational demanding approach of LCMs to assess the advantages of high-resolution cloud models and further to estimate errors of traditional parameterizations.publishedVersio

Can Areawide Building Retrofitting Affect the Urban Microclimate? An LES Study for Berlin, Germany

In this work, we investigate the effect of areawide building retrofitting on summertime, street-level outdoor temperatures in an urban district in Berlin, Germany. We perform two building-resolving, weeklong large-eddy simulations: one with nonretrofitted buildings and the other with retrofitted buildings in the entire domain to meet today’s energy efficiency standards. The comparison of the two simulations reveals that the mean outdoor temperatures are higher with retrofitted buildings during daytime conditions. This behavior is caused by the much smaller inertia of the outermost roof/wall layer in the retrofitting case, which is thermally decoupled from the inner roof/wall layers by an insulation layer. As a result, the outermost layer heats up more rigorously during the daytime, leading to increased sensible heat fluxes into the atmosphere. During the nighttime, the outermost layer’s temperature drops down faster, resulting in cooling of the atmosphere. However, as the simulation progresses, the cooling effect becomes smaller and the warming effect becomes larger. After 1 week, we find the mean temperatures to be 4 K higher during the daytime while the cooling effects become negligible.publishedVersio

An Improved Surface Boundary Condition for Large-Eddy Simulations Based on Monin–Obukhov Similarity Theory: Evaluation and Consequences for Grid Convergence in Neutral and Stable Conditions

Monin–Obukhov similarity theory is used in large-eddy simulation (LES) models as a surface boundary condition to predict the surface shear stress and scalar fluxes based on the gradients between the surface and the first grid level above the surface. We outline deficiencies of this methodology, such as the systematical underestimation of the surface shear stress, and propose a modified boundary condition to correct for this issue. The proposed boundary condition is applied to a set of LES for both neutral and stable boundary layers with successively decreasing grid spacing. The results indicate that the proposed boundary condition reliably corrects the surface shear stress and the sensible heat flux, and improves grid convergence of these quantities. The LES data indicate improved grid convergence for the surface shear stress, more so than for the surface heat flux. This is either due to a limited performance of the Monin–Obukhov similarity functions or due to problems in the LES model in representing stable conditions. Furthermore, we find that the correction achieved using the proposed boundary condition does not lead to improved grid convergence of the wind-speed and temperature profiles. From this we conclude that the sensitivity of the wind-speed and temperature profiles in the LES model to the grid spacing is more likely related to under-resolved near-surface gradients and turbulent mixing at the boundary-layer top, to the SGS model formulation, and/or to numerical issues, and not to deficiencies due to the use of improper surface boundary conditions.publishedVersio

Modeling of land–surface interactions in the PALM model system 6.0: land surface model description, first evaluation, and sensitivity to model parameters

In this paper the land surface model embedded in the PALM model system is described and evaluated against in situ measurements at Cabauw, Netherlands. A total of 2 consecutive clear-sky days are simulated, and the components of surface energy balance, potential temperature, humidity, and horizontal wind speed are compared to observations. For the simulated period, components of the energy balance are consistent with daytime and nighttime observations, and the daytime Bowen ratio also agrees fairly well with observations. The model simulates a more stably stratified nocturnal boundary layer than the observations, and near-surface potential temperature and humidity agree fairly well during the day. Moreover, a sensitivity analysis is performed to investigate dependence of the model on land surface and soil specifications, as well as atmospheric initial conditions, because they represent a major source of uncertainty in the simulation setup. It is found that an inaccurate estimation of leaf area index, albedo, or initial humidity causes a significant misrepresentation of the daytime turbulent sensible and latent heat fluxes. During the night, the boundary-layer characteristics are primarily affected by surface roughness and the applied radiation schemes.publishedVersio

Evaluation of the dynamic core of the PALM model system 6.0 in a neutrally stratified urban environment: comparison between LES and wind-tunnel experiments

We demonstrate the capability of the PALM model system version 6.0 to simulate neutrally stratified urban boundary layers. Our simulation uses the real-world building configuration of the HafenCity area in Hamburg, Germany. Using PALM's virtual measurement module, we compare simulation results to wind-tunnel measurements of a downscaled replica of the study area. Wind-tunnel measurements of mean wind speed agree within 5 % on average while the wind direction deviates by approximately 4∘. Turbulence statistics similarly agree. However, larger differences between measurements and simulation arise in the vicinity of surfaces where building geometry is insufficiently resolved. We discuss how to minimize these differences by improving the grid layout and give tips for setup preparation. Also, we discuss how existing and upcoming features of PALM like the grid nesting and immersed boundary condition help improve the simulation results.publishedVersio

The Departure from Mixed-Layer Similarity During the Afternoon Decay of Turbulence in the Free-Convective Boundary Layer: Results from Large-Eddy Simulations

This study analyses the departure of the velocity-variances profiles from their quasi-steady state described by the mixed-layer similarity, using large-eddy simulations with different prescribed shapes and time scales of the surface kinematic heat flux decay. Within the descriptive frames where the time is tracked solely by the forcing time scale (either constant or time-dependent) describing the surface heat flux decay, we find that the normalized velocityvariances profiles from different runs do not collapse while they depart from mixed-layer similarity. As the mixed-layer similarity relies on the assumption that the free-convective boundary layer is in a quasi-equilibrium, we consider the ratios of the forcing time scales to the convective eddy-turnover time scale. We find that the normalized velocity-variances profiles collapse in the only case where the ratio (r) of the time-dependent forcing time scale to the convective eddy-turnover time scale is used for tracking the time, supporting the independence of the departure from the characteristics of the surface heat flux decay. As a consequence of this result, the knowledge of r is sufficient to predict the departure of the velocity variances from their quasi-steady state, irrespective of the shape of the surface heat flux decay. This study highlights the importance of considering both the time-dependent forcing time scale and the convective eddy-turnover time scale for evaluating the response of the free-convective boundary layer to the surface heat flux decay

Importance of radiative transfer processes in urban climate models: A study based on the PALM 6.0 model system

Including radiative transfer processes within the urban canopy layer into microscale urban climate models (UCMs) is essential to obtain realistic model results. These processes include the interaction of buildings and vegetation with shortwave and longwave radiation, thermal emission, and radiation reflections. They contribute differently to the radiation budget of urban surfaces. Each process requires different computational resources and physical data for the urban elements. This study investigates how much detail modellers should include to parameterize radiative transfer in microscale building-resolving UCMs. To that end, we introduce a stepwise parameterization method to the Parallelized Large-eddy Simulation Model (PALM) system 6.0 to quantify individually the effects of the main radiative transfer processes on the radiation budget and on the flow field. We quantify numerical simulations of both simple and realistic urban configurations to identify the major and the minor effects of radiative transfer processes on the radiation budget. The study shows that processes such as surface and vegetation interaction with shortwave and longwave radiation will have major effects, while a process such as multiple reflections will have minor effects. The study also shows that radiative transfer processes within the canopy layer implicitly affect the incoming radiation since the radiative transfer model is coupled to the radiation model. The flow field changes considerably in response to the radiative transfer processes included in the model. The study identified those processes which are essentially needed to assure acceptable quality of the flow field. These processes are receiving radiation from atmosphere based on the sky-view factors, interaction of urban vegetation with radiation, radiative transfer among urban surfaces, and considering at least single reflection of radiation. Omitting any of these processes may lead to high uncertainties in the model results.publishedVersio

A nested multi-scale system implemented in the large-eddy simulation model PALM model system 6.0

Large-eddy simulation (LES) provides a physically sound approach to study complex turbulent processes within the atmospheric boundary layer including urban boundary layer flows. However, such flow problems often involve a large separation of turbulent scales, requiring a large computational domain and very high grid resolution near the surface features, leading to prohibitive computational costs. To overcome this problem, an online LES–LES nesting scheme is implemented into the PALM model system 6.0. The hereby documented and evaluated nesting method is capable of supporting multiple child domains, which can be nested within their parent domain either in a parallel or recursively cascading configuration. The nesting system is evaluated by first simulating a purely convective boundary layer flow system and then three different neutrally stratified flow scenarios with increasing order of topographic complexity. The results of the nested runs are compared with corresponding non-nested high- and low-resolution results. The results reveal that the solution accuracy within the high-resolution nest domain is clearly improved as the solutions approach the non-nested high-resolution reference results. In obstacle-resolving LES, the two-way coupling becomes problematic as anterpolation introduces a regional discrepancy within the obstacle canopy of the parent domain. This is remedied by introducing canopy-restricted anterpolation where the operation is only performed above the obstacle canopy. The test simulations make evident that this approach is the most suitable coupling strategy for obstacle-resolving LES. The performed simulations testify that nesting can reduce the CPU time up to 80 % compared to the fine-resolution reference runs, while the computational overhead from the nesting operations remained below 16 % for the two-way coupling approach and significantly less for the one-way alternative.publishedVersio